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Abstract 

We consider an anisotropic homogenized composite medium (HCM) arising from isotropic par- 
ticulate component phases based on elhpsoidal geometries. For cubically nonlinear component 
phases, the corresponding zeroth-order strong-permittivity-fluctuation theory (SPFT) (which is 
equivalent to the Bruggeman homogenization formalism) and second-order SPFT are established 
and used to estimate the constitutive properties of the HCM. The relationship between the com- 
ponent phase particulate geometry and the HCM constitutive properties is explored. Significant 
differences are highlighted between the estimates of the Bruggeman homogenization formalism and 
the second-order SPFT estimates. The prospects for nonlinearity enhancement are investigated. 



1 Introduction 

The constitutive properties of a homogenized composite medium (HCM) are determined by both 
the constitutive properties and the topological properties of its component phases [l]-[5]. In partic- 
ular, component phases based on nonspherical particulate geometries may give rise to anisotropic 
HCMs, despite the component phases themselves being isotropic with respect to their electromag- 
netic properties. Such geometrically-derived anisotropy has been extensively characterized for linear 
dielectric HCMs [6]-[8] and more general bianisotropic HCMs [8]-[10]. For weakly nonlinear HCMs, 
the role of the component phase particulate geometry was emphasized recently in this journal by 
Goncharenko, Popelnukh and Venger [11], using an approach founded on the mean-field approxi- 
mation. However, their analysis was restricted to the Maxwell Garnett homogenization formalism 
[3, 12]. A more comprehensive study is communicated here based on the strong-permittivity- 
fluctuation theory (SPFT) [13]. In contrast to the aforementioned Maxwell Garnett approach [11], 
the SPFT approach (i) incorporates higher-order statistics to describe the component phase dis- 
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tributions; (ii) is not restricted to only dilute composites; and (iii) is not restricted to only weakl^ 
nonspherical particulate geometries. 

The early development of the SPFT concerned wave propagation in continuous random medi- 
ums [14, 15], but more recently the theory has been applied to the estimation of HCM constitutive 
parameters [16, 17]. The SPFT represents a significant advance over conventional homogcnization 
formalisms, such as the Maxwell Garnett approach and the Bruggeman approach [3, 6], through 
incorporating a comprehensive description of the distributional statistics of the HCM component 
phases. In estimating the constitutive parameters of an HCM, the SPFT employs a Feynman- 
diagrammatic technique to calculate iterative refinements to the constitutive parameters of a com- 
parison medium; successive iterates incorporate successively higher-order spatial correlation func- 
tions. It transpires that the SPFT comparison medium is equivalent to the effective medium of 
the (symmetric) Bruggeman homogenization theory [20, 21]. In principle, correlation functions of 
arbitrarily high order can be accommodated in the SPFT. However, the theory is most widely- 
implemented at the level of the bilocal approximation (i.e., second-order approximation), wherein 
a two-point covariance function and its associated correlation length L characterize the component 
phase distributions. As indicated in figure 1, coherent interactions between pairs of scattering 
centres within a region of linear dimensions L arc incorporated in the bilocal SPFT; scattering 
centres separated by distances much greater than L are assumed to act independently. Thereby, 
the SPFT provides an estimation of coherent scattering losses, unlike the Maxwell Garnett and 
Bruggeman homogenization formalisms. Notice that the bilocally-approximated SPFT gives rise 
to the Bruggeman homogenization formalism in the limit L — > [21]. 

The SPFT has been widely applied to linear homogenization scenarios, where generalizations^ 
have been developed for anisotropic dielectric [18, 19], isotropic chiral [20] and bianisotropic [21, 
22] HCMs. Investigations of the trilocally-approximated SPFT for isotropic HCMs have recently 
confirmed the convergence of the second-order theory [17, 23, 24]. In the weakly nonlinear regime, 
developments of the bilocally-approximated SPFT have been restricted to isotropic HCMs, based 
on spherical component phase geometry [16, 17, 24]. The present study advances the nonlinear 
SPFT through developing the theory for cubically nonlinear, anisotropic HCMs. Furthermore, it 
is assumed that the component phases are composed of electrically-small ellipsoidal particles. The 
relationship between the HCM constitutive parameters and the underlying particulate geometry of 
the component phases is investigated via a representative numerical example. 

In our notational convention, dyadics are double underlined whereas vectors are in bold face. 
The inverse, adjoint, determinant and trace of a dyadic A are denoted by , det 



A 



and 



tr 



A 



respectively. The identity dyadic is represented by J_ . The ensemble average of a quantity 
ijj is written as ( ^ ). The permittivity and permeability of free space (i.e., vacuum) are given by eo 
and fio, respectively; ko = LOy/e^JI^ is the free-space wavenumber while uj is the angular frequency. 

^The generalized SPFT is referred to as the strong-property-fluctuation theory. 



2 Homogenization generalities 
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2.1 Component phases 

Consider the homogenization of a two-phase composite with component phases labehed as a and 
b. The component phases are taken to be isotropic dielectric mediums with permittivities 

e£ = e« + x^|E£p, {i = a,b), (1) 

where e^o is the hnear permittivity, xe is the nonlinear susceptibility, and | ^ is the electric field 
developed inside a region of phase i by illumination of the composite medium. We assume weak 
nonlinear ity; i.e., | e^o I S> | | | p. Notice that such electrostrictive mediums as characterized 
by (1) can induce Brillouin scattering which is often a strong process [25]. The component phases 
a and b are taken to be randomly distributed as identically-orientated, conformal ellipsoids. The 
shape dyadic 

a = \ diag(C/„ Uy, U,), Uy, U, > 0), (2) 

^UxUyUz 

parameterizes the conformal ellipsoidal surfaces as 

re{eA) = vE'r{e,(j)), (3) 

where f {0, cp) is the radial unit vector specified by the spherical polar coordinates 9 and (p. Thus, 
a wide range of ellipsoidal particulate shapes, including highly elongated forms, can be accommo- 
dated. The linear ellipsoidal dimensions, as determined by r), are assumed to be sufficiently small 

that the electromagnetic long-wavclcngth regime pertains. 

In the SPFT, statistical moments of the characteristic functions 



1, r G Ve, 

{i = a,b), (4) 

0, r ^ Ve, 



are utilized to take account of the component phase distributions. The volume fraction of phase £, 
namely fi , is given by the first statistical moment of ; i.e., ( $^(r) ) = fe ■ Clearly, fa + fb = 
1. The second statistical moment of provides a two-point covariance function; we adopt the 
physically-motivated form [26] 

{^e{r)<^e{r'))=fe [1 + ( /, - 1 ) H ( a - L ) ] , (5) 

where 7i is the Heavisidc function (i.e., 7i{x) = ^l^^iu) dy where 5 is the Dirac delta function), 
(7 = I C/~^ • R| with R = r — r', and L > is the correlation length. The specific nature of the 
covariance function has been found to exert little influence on the SPFT estimates for linear [19] 
and weakly nonlinear [17] HCMs. 

2.2 Homogenized composite medium 

Let 'FiHCM denote the spatially-averaged electric field in the HCM. In this communication we derive 
the estimate 

= IfcaO+^fejEz/CMl' (6) 

= diag (ef„o, eLo' eLo) + diag (x^, xL' xla) I ^hcm \^ (7) 



of the HCM permittivity. The bilocally- approximated SPFT is utilized (hence the subscripts ba i^ 
(6), (7) ). Note that the Bruggeman estimate of the HCM permittivity, namely 

^Br = ^BrO+KsJ^HCM^ (8) 

= diag (e|^o > ^BrO ' ^BrO ) + diag {xbv > Xbt ^ Xbv) I E hcm f, (9) 

characterizes the comparison medium which is adopted in the bilocally-approximated SPFT [21]. 
As the Bruggeman homogenization formalism — in which the component phases a and b are treated 

symmetrically [6] — provides the comparison medium, the SPFT homogenization approach (like 
the Bruggeman formalism) is applicable for all volume fractions /„ G (0, 1). 



2.3 Depolairization and polarizability dyadics 

The depolarization dyadic ^ is a key element in both Bruggeman and SPFT homogenizations. It 
provides the electromagnetic response of a ^-shaped exclusion volume, immersed in a homogeneous 
background, in the limit rj ^ 0. For the component phases described by (1) and (2), we find [27, 28] 



D = d(t> r de suiO \ r ^ tA\ 

— luiir Jo Jo yti gj^^ 'A ~ I 



wherein 



A = diag 



sin^ 6 cos^ ^ siv? 9 sin^ 4> cos^ 9 



Uy 



(10) 



(11) 



The integrations of (10) reduce to elliptic function representations [29]. In the case of spheroidal 

particulate geometries, hyperbolic functions provide an evaluation of _D [27], while for the degenerate 
isotropic case Ux = Uy = Uz we have the well-known result Q, = (l/iwS) [30]. We express Q_ as 
the sum of linear and weakly nonlinear parts 



D^ = D^o+^i\'^HCM 



(12) 



with 



^0 



^1 



f2lT 




/ d( 


p d9 sine 1 


/o 





tr 



1 



2tt 



d9 sin 9 



tr 



X ' A 

= Br 



I V 



(13) 



(14) 



A convenient construction in homogenization formalisms is the polarizability dyadic X ^, defined 



as 



-lUO 



Br J =i 



{e = a,b), 



where 



I_ + iuj D' ( Q X — e 



Let us proceed to calculate the linear and nonlinear contributions in the decomposition 

'2 {i = a,b). 



(15) 
(16) 

(17) 



Under the assumption of weak nonlinear ity, we express (16) in the form 
with Hnear term 
and nonlinear term 
The local field factor 
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(18) 
(19) 



{aexa-K^J+Ri' i^eoL-g^,,)] • (20) 



d\Bf 



d E 



HCM 



{i = a,b), 



(21) 



has been incorporated in deriving (18)-(20), via the Maclaurin series expansion = e^Q+g^ Xi I ^HCM 
An appropriate estimation of the local field factor is provided by [31] 



tr 



Thus, the inverse of £^ is given as 



(22) 



(23) 



wherein 

A, 

and 



det 



r 

=£0 



diag (^r^^rfo + r^o^li' ^firfo + r|orfi, r^^rf o + r^o^fi ) ~ peL 



(24) 



Pi — ^eo^^o^ei + ^eo^n^eo + ^ei^^o^eo ■ (25) 
Combining (23) and (24) with (15), and separating linear and nonlinear terms, provides 



X 



eo 



iuj Qol - e, 



BrO J =eO 



-ILO 



Lei + {^eoL-g 



BrO J =£} 



{i = a,b). (26) 



2.4 Bruggeman homogenization 

The Bruggeman estimates of the HCM linear permittivity e^^^ and nonlinear susceptibility x ^ 
are delivered through solving the nonlinear equations [3, 6, 31] 



(j = o,i). 



(27) 



Recursive procedures for this purpose provide the p*^ iterates [4, 24] ^ 

(28) 



in terms of the (p — 1)*'^ iterates, wherein the operators are defined by 

{ i baa ] = [f'^ + E^o ) • ( £ ad + £ 60 ) ' 

^^{^ba} = {/a[5aXaE;;+(^«oi-|6ao) * 4 J + A [^6 X6 I^o + (^"0 I - | J '46]} 

• (/aE;d+/6E;oT'' (29) 

while suitable initial values are given by 

= ( f„ f^n A- fu fur, ^ 7^ I 

(30) 



ifeaO = {fa€aO + fb^bo) L 
K^J^^ = ifaXa + fbXb) L 



3 The bilocally-approximated SPFT 

The bilocally-approximated SPFT estimate of the HCM permittivity dyadic, as derived elsewhere 
[21], is given by 

the mass operator term 

^ba=Ua-^b) Ua-^b) (^2) 

is specified in terms of the principal value integral 

m = r [ ^B,(R), (33) 

Ja<L 

with G^^(R) being the unbounded dyadic Green function of the comparison medium. Here we 
develop expressions for the linear and nonlinear contributions oi e^^, appropriate to the component 

phases specified in §2. 

Under the assumption of weak nonlinearity, we express W_ = W_q + W_^ \ ^hcm P; integral 
expressions for and are provided in the Appendix. Thereby, the linear and nonlinear 
terms in the mass operator decomposition S = S^^^ + S^^^^ | Ehcm P are given as 

i6aO=(^aO-^6o) ' ^0 ' (^aO " ^60 ) ' (3^) 
i6al=2(^„0-^6o) -Eo- Ual-^6l) + (^aO-^6o) * El* U„0 " ^60 ) ' (^5) 

respectively, correct to the second order in | Ehcm \ ■ Now, let us introduce the dyadic quantity 

= I + i6a '^ = 00+01 lE/fCM I', (36) 



such that 
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We may then express the inverse dyadic in the form 



with nonhnear part 



n 



det 



where 



(37) 
(38) 



(39) 



diag (offjg + n^ni, ofog + n^o^f, n^n^ + n^nf^ - lyg'/A , (40) 



(41) 



Thus, the linear and nonhnear contributions of the SPFT estimate e^^^ are dehvered, respectively, 
as 



-baO 



— • E 



iBrO i^^=0 



baO 



4a ^Br 



^{^''^bal+U'Mbao) 



(42) 
(43) 



4 Numerical results and discussion 



Let us explore the HCM constitutive parameter space by means of a representative numerical 
example: Consider the homogenization of a cubically nonlinear phase a with linear permittivity 
Sao = 2eo and nonlinear susceptibility Xa = 9.07571 x IQ-^^eo m^V"^ (= 6.5 x 10~^ esu) and a linear 
phase b with permittivity ei, = e;,o = 12eo. Note that the selected nonlinear susceptibility value 
corresponds to that of gallium arsenide [25], while selected the linear permittivity values are typical 
of a wide range of insulating crystals [32]. We assume the elhpsoidal component phase topology 
specified by Ux = I, Uy = 3 and Uz G [0.5, 15]. The angular frequency u is fixed at 27r x 10-'^'^ rad s-"*^ 
for all calculations reported here. 

The Bruggeman estimates of the HCM relative linear and nonlinear constitutive parameters 
are plotted in figure 1 as functions of /„ and Uz- The calculated constitutive parameters presented 
in figure 1 are consistent with those calculated by Lakhtakia and Lakhtakia [31] in a study per- 
taining to the Bruggeman homogenization of ellipsoidal inclusions with a host medium comprising 
spherical particles. The linear parameters follow an approximately linear progression between their 
constraining values at = and /„ = 1. Furthermore, for the range Uz G [0.5,15], the linear 
parameters are largely (but not completely) independent of the particulate geometry of the com- 
ponent phases. This is in contrast to the nonlinear parameters which are acutely sensitive to Uz- 
Of special significance is the nonlinearity enhancement (i.e., the manifestation of a higher degree 
of nonlinear susceptibility in the HCM than is present in its component phases) which is partic- 
ularly observed at high values of Uz for xhr a^*^ values of Uz for Xbv This phenomenon 
and its possible technological exploitation are described elsewhere [3, 16, 17, 24, 31, 33]. In order 



to best consider nonlinearity enhancement, we fix the shape parameter Uz = 15 for all remaining 

calculations. 

We turn our attention now to the bilocally-approximated SPFT calculations. Let 



,nr^ = 'baO ^BrO ^ ^n. ^ Xfc„ Xsr ^ {n = X,y,z). (44) 

Co Xa 

The SPFT estimates of the HCM relative linear constitutive parameters e^^^^'^' and nonlinear con- 
stitutive parameters xta''^^'^^ plotted in figures 2 and 3, respectively, as functions of /„ and koL. 
Significant differences are clear between the Bruggeman-estimated values and the SPFT-estimated 
values: The SPFT estimates of linear constitutive parameters provide an additive correction to the 
corresponding Bruggeman parameters, whereas for the nonlinear constitutive parameters the SPFT 
estimates provide a subtractive correction to the corresponding Bruggeman parameters. Further- 
more, the magnitudes of these differences exhibit local maxima which occur at progressively higher 
values of /„ as one compares the constitutive parameter components aligned with the x, y and z 
coordinate axes, respectively. This trend holds for both the real and the imaginary parts of both 
the linear permittivity and the nonlinear susceptibility parameters. However, it is less pronounced 
for the nonlinear constitutive parameters. 

Coherent interactions between scattering centres enclosed within a region of linear dimensions 
L are accommodated in the bilocally-approximated SPFT via the two-point covariance function 
(5) (see figure 1). Thus, since neither component phase a nor component phase h is dissipative, the 
nonzero imaginary parts of the SPFT constititutive parameters in figures 2 and 3 are attributable 
entirely to scattering losses. Furthermore, the magnitudes of the imaginary parts of the constitutive 
parameters are observed in figures 2 and 3 to increase as L increases, due to the actions of greater 
numbers of scattering centres becoming correlated. 



5 Concluding remarks 

The bilocally-approximated SPFT for weakly nonlinear isotropic HCMs, based on spherical partic- 
ulate geometry, has been recently established [16, 17, 24]. In the present study we further advance 
the theory through considering anisotropic, cubically nonlinear HCMs, arising from isotropic com- 
ponent phases with ellipsoidal particulate geometries. Significant differences between the bilocally- 
approximated SPFT (i.e., second-order theory) and the Bruggeman homogenization formalism (i.e., 
zeroth-order theory) — which depend upon the underlying particulate geometry — have emerged. 
In particular, nonlinearity enhancement is predicted to a lesser degree with the SPFT than with the 
Bruggeman homogenization formalism. The importance of taking into account the distributional 
statistics of the HCM component phases is thereby further emphasized. 
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Consider the principal value integral term W_ (33) which was expressed as the sum W ^ + 
W ^ I 'Ehcm P in §3. Here we develop expressions for the linear component and the nonlinear 
component W-i, appropriate to the homogenization scenario of §2. 

Let us begin with the following straightforward specialization of the evaluation of W_ for bian- 
isotropic HCMs [21, 22] ~ 



W = 

= 27r2 • 



ILU 



{q/uftc + {qluftB + tA 



sin qL 



— L cos qL 



(45) 



For the weakly nonlinear homogenization outlined in §2, the scalar terms tA, ts and tc in (45), 
along with their linear and nonlinear decompositions, are given by 



I Br 



tAO + tAl I E HCM I , 



tr 



adj 
= Br 



tr 



A 



I = tso + tBl I ^HCM 



tc = Ho tr 



A 



tr 



tco + tci I E HCM 



(46) 

(47) 
(48) 



wherein 



tAO = ^0 det 



* Al — ^0 ( XBr <^iJrO *'BrO + ^BrO Xsr ^BrO + ^BrO ^SrO XSr ) ' (^9) 



tr 
T 



tr 



tBi=nl {trfX'i] - (ti" 
X = diag ( XBr ^BrO + ^BrO XSr) Xfir ^SrO + ^BrO XSr i XSr ^SrO + ^BrO XBr ) ' 



tr 



4])}. 



ico = A«o tr 



A 



tr 



tci = Mo tr 



tr 



=±Jr — 



(50) 

(51) 
(52) 
(53) 



Similarly, the dyadic quantities a and P in (45), along with their linear and nonlinear decomposi- 
tions, are given by 



a = Ho 



(2|B,-tr[|5j l) -i-tr 



^Br'A 



tB. 
- /Xo — tr 

tc 



A 



£ = A^olB;-Mo^tr 



A 



4 = £o+£JEijcMr> 



(54) 
(55) 



with 



a. 



a. 



/3 

= 



Mo 



Mo 



(2^ 



Br- 



tr 



zBrO 



^Br 



X ' ^ 

= Br = 



- ^0 T^tr \a \ a, 
tco '- 

tBl — {t BO /tco) tci . 

Mo 7 tr 

tco 



A 



3 adj 
M-o EBrO 



tAO . 
Ho- — tr 
ceo 



A 



A. 



a 3 tAl — {tAO tco) tci . 

=MoX-Mo tr 

= 1 ~ tco 



A 



A. 



(56) 

(57) 
(58) 



In the long- wavelength regime, i.e., | d± | ^ 1 , the application of residue calculus to (45) deliver 



W 



Am 



Jo Jo 



de 



sin( 



3A a;2 



a + i — t 



where we have introduced 



A = sJtl-AtAtc = Ao + Ai I ^HCM |^ 
2 -is ± A 



K±= UJ 



2tc 



K0± + Kl± ^HCM , 



d± = L^/k± = do± + di± I E HCM 
with linear and nonlinear parts 



An 



2 ~tBO ± Ao 



Ai 



'^"^^^ 2 too 



tBotsi — 2 {tAltco + tAptci ) 

Ao 

2 ( -tBi ± Ai ) - 2tci (ko± V) 



2ico 



The linear and nonlinear components of W_ are thereby given as 



W = 

= 



27r 



and 



{ 1 



fafb(^ 

Airi Jo 



dO sm9 X 



dO 



3Ao 



/'27r /"TT 

/ / 

JO Jo 

r 3 r 

a , + rg /3 H — ^ do+<^i+ (1 + «<^o+) — <io_cZi_ (1 + zc/o-) 
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The integrals (66) and (67) are straightforwardly evaluated by standard (e.g., Gaussian) nu- 
merical methods [34]. In the degenerate isotropic case = Uy = Uz, the integrals (66) and (67) 
yield the analytic results of [17]. 



Figure 1: Schematic diagram illustrating the bilocally-approximated SPFT for ellipsoidal com- 
ponent phase geometry: pair-wise scattering interactions are accommodated between ellipsoidal 
scattering centres contained within an ellipsoidal correlation region of linear dimensions L. 
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Figure 2: HCM relative linear permittivity and nonlinear susceptibility parameters calculated using 
the Bruggeman homogenization formalism. Component phase parameter values: Cao = 2eo, Xa = 
9.07571 X 10-^2g^jj^2Y-2^ ^ ^ ;L2eo, C/^ = 1 and C/^ = 3. 




Figure 3: Real and imaginary parts of the HCM linear permittivity parameters calculated using the 
SPFT homogenization formalism. Component phase parameter values: eao = 2eo, Xa = 9.07571 x 
10-^2g^^2Y-2^ ^ ^ ^2eo, Ua; = l,Uy = 3 and = 15. 
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Figure 4: As figure 3 but for the HCM nonlinear susceptibility parameters. 



